Anisotropic conductivity tensor on a half-filled high Landau level 



O 

o 

(N 

>>: 



oo 



I.S. Burmistrov 

Landau Institute for Theoretical Physics, Kosygina str. 2, 117940 Moscow, Russia and 
Institute for Theoretical Physics, University of Amsterdam, 
Valckenierstraat 65, 1018XE Amsterdam, The Netherlands 

We study two-dimensional interacting electrons in a weak perpendicular magnetic field with the 
filling factor v 3> 1 and in the presence of a quenched disorder. As it is known, the unidirectional 
charge density wave state can exist near a half-filled high Landau level at low temperatures if 
disorder is weak enough. We show that the existence of the unidirectional charge density wave state 
at temperature T < T c where T c is the transition temperature leads to the anisotropic conductivity 
tensor. We find that the anisotropic part of conductivity tensor is proportional to (T c — T)/T c below 
the transition in accordance with the experimental findings. The order parameter fluctuations wash 
out the mean-field cusp at T = T c and the conductivity tensor becomes anisotropic even above the 
mean-field transition temperature T c . 
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PACS numbers: 72.10.-d, 73.20.Dx, 73.40.Hm 

I. INTRODUCTION 

A two-dimensional electron gas (2DEG) in a perpen- 
dicular magnetic field (H) have remained a subject of 
intensive studies, both theoretical and experimental, for 
several decades^. Recently the phenomenon of sharp 
anisotropy of magnetoresistance near half-filled high Lan- 
dau levels with Landau level index N = 2,3,4 at low 
temperatures have been discovered^. Right away, the 
magnetoresistance anisotropy have been related with the 
possible existence of unidirectional charge density wave 
(UCDW) state near half-filling of a high Landau level 
that theoretical results of RefsA^ had been predicted 
in the limit N 3> 1. Later a bunch of theoretical and ex- 
perimental papers have appeared^. However, only a few 
theoretical researches have been concerned the transport 
properties of the system whereas the only experimentally 
measured quantity is magnetoresistance. 

In the limit of zero temperature the unidirectional 
charge density wave has strongly non-sinusoidal profile 
with well-defined edges^ and the conductivity tensor can 
be evaluated in the framework of an effective theory for 
edge excitations^. The result of Ref£ is in a qualitative 
agreement with measured magnetoresistance behavior in 
the limit of zero temperature^. Later the "semicircle 
law" of Dykhnc-RuziniS have also been extended to the 
case of the anisotropic conductivity tensor—. 

At present a thorough analysis of temperature (T) de- 
pendence of magnetoresistance near the transition tem- 
perature T c from the liquid state to the UCDW state is 
absent. The main objective of the paper is to present 
such analysis for the domain \T C — T\/T c -C 1 where the 
expansion in the UCDW order parameter is legitimate. 
The effect of the order parameter fluctuations that are 
enhanced near the phase transition on the conductivity 
tensor both below and above the transition temperature 
T c is also investigated. 

We assume the presence of a weak random potential 
created by impurities near 2DEG such that the Landau 
level broadening 1 /r is much less the spacing u>h between 



Landau levels, 1/t <C u>jj. Here the u>h is cyclotron 
frequency u>h = eH/m with e and m being the electron 
charge and band mass respectively (we use the units with 
H = 1, c = 1, and ks = 1). 

One among the main results of the paper is the fact 
that the conductivity tensor a a b (we measure conduc- 
tivity in units of e 2 /h) of two-dimensional electrons in 
the UCDW state acquires anisotropic part cr^ nis ' ) propor- 
tional to temperature deviation from T c , i.e. cr^ nls ^ oc 
(T c — T)/T c . Other main result is that in the vicin- 
ity of the transition temperature T c there is additional 
anisotropic contribution <5cr^ ms_f ^ cx (\T C — T\/T c )~ 3 / 2 
that we will refer as fluctuational to the conductivity ten- 
sor a^J of the liquid state due to UCDW order parameter 
fluctuations. 

We start out with an introduction to the formalism 
that mainly follows one introduced in the previous pa- 
per—. In section [H] the effective action of "three level 
model" is developed. The conductivity tensor is evalu- 
ated in section 11111 The effect of order parameter fluc- 
tuations is investigated in section IIVI In section the 
results obtained are discussed in relation with recent ex- 
periments^. 

Some of the results of the present paper have been 
published in a brief form in Ref— . 



II. "THREE LEVEL" MODEL 

A. Introduction 

To start out consider the system of two-dimensional 
interacting electrons in the presence of a random poten- 
tial Vdis(r) and perpendicular magnetic field H. The pa- 
rameter that characterizes the strength of the Coulomb 
interaction is r s = \[2e 2 jevp with vp being the Fermi 
velocity and e the dielectric constant of a media. We as- 
sume that the Coulomb interaction between the electrons 
is weak, r s <l, and the magnetic field obeys the condi- 
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tion Nr s ^> 1 where the Landau level index N = \v/2\ is 
the integer part of half of filling factor v. In addition, we 
assume that the Landau level broadening 1 /r is not too 
small, 1/r > luhN^ 1 \n^/2r s N. In this case it is possible 
to construct an effective field theory for electrons on the 
highest partially filled Landau level by integrating out all 
other degrees of freedomi^i^. 

Also, we consider the case when the Landau lev- 
els are spin-resolved that occurs according to Refii^ if 
1/r < A cx = (r s u>H/nV2) ln2\/2/r s . Therefore, the 
Landau level broadening should be restricted from below 
and from above as 

^lnV2r s iV«I«^ln^. (1) 

Throughout the paper we characterize the UCDW 
state by the order parameter A(Q). All calculations are 
performed under the assumption Nr 2 . 3> 1. In this case 
the Hartree-Fock approximation is well justified^ more- 
over it is known that the corrections to it are small in the 
parameter a B /l H = 1/Ar 2 <C 1, where a B = e/me 2 is 
the Bohr radius and Ih = 1/ the magnetic length. 

B. Formalism 



magnetic field with vector-potential a, 

H= -'-HV-eA-ea) 2 . (6) 
2m 

As usual, we assume the white-noise distribution for 
the random potential, 

V[V dis (r)} = -±= exp (~ J dvVl(r)\ . (7) 

In order to average over disorder we introduce N r repli- 
cated copies of the system^ labelled by the replica in- 
dices a = 1, • • ■ ,N r . 

It is convenient to rewrite one-particle hamiltonian © 
with a help of covariant derivative 

D = V - ieA, (8) 

in order to extract the time-dependent vector potential 
a(V„), 

W = -^D 2 +A>„), (9) 

9 

K{v n ) = a(f„)D + — a(i^„_ m ) a(i/ m ). 

m. /.m — ' 



Thermodynamic potential of the system in hand can 
be written as 

Sl = -^r fvfoifl J V[V dis ]V[V dis ] expS^V.^dis]. 

\ (2) 
where action S[tp, if>, Vdi s ] is written in Matsubara repre- 
sentation 

E /rfrdr'Cf(r)C;%(r) 
xU Q (r-r>)^(r)i>Z + (r'). (3) 



Here, ipZ^{ r ) an( i ^"^(r) are annihilation and creation 
electron operators. T stands for temperature, fj. chemi- 
cal potential, a and a spin indices, Lu n = irT(2n + 1) 
fermionic frequency whereas v n = 2i:Tn bosonic one. 
Matrix H is defined as 

w = J>(i/ n )C (4) 

CtV n 

with matrices 

{Infki = S aP 5^5 k -i, n (5) 

being U(l) generators. One-particle hamiltonian H. de- 
scribes a two-dimensional electron in constant perpendic- 
ular magnetic field H = e a bd a Ab and in time-dependent 



C. Effective action of "three-level" model 

To investigate the thermodynamic properties of elec- 
trons on the Nth Landau level one can integrate out 
electrons on all other Landau levels^. However, to find 
conductivity tensor projection on the single Ath Landau 
level is not appropriate because of covariant derivative 
D has non-zero matrix elements only for transitions be- 
tween adjacent Landau levels. It is necessary therefore 
to consider not only the iVth Landau level alone but two 
adjacent ones, the (A— l)th and (iV+l)th Landau levels. 

Extending the projection to the Ath Landau level only 
of RefsiMii^, we obtain effective action for electrons on 
the (A-l)th, Ath, and (A + l)th Landau level as follows 

xCrW-f E /drdr'Cf(r)C;%,(r) 

x[/ scr (r-r')C^(r)C;:^(r')- (10) 



Here i/ , 2'°'( r ) ano - V'S,',^ 1- ) arc annihilation and creation 
operators of an electron on the (A — l)th, Ath, and (A + 
l)th Landau levels, 

N+l N+l 

c»= E c«. Wv= E «(>•)• 

p=N-l p=N-l 

(11) 
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The screened electron-electron interaction U SCI (r) has the 
following Fourier transform 



£4cr(<?) 



2ne 2 
eq 



qa B \ 6uj h t 



(12) 



[1- J Q 2 (qR c )-2J?(qR c )) 



It is different from one obtained in Refi^ii. The reason 
for that is exclusion of contributions from the (N — l)th 
and (N + l)th Landau level from the polarization oper- 
ator. 

Effective action <|10|) was obtained under assumptions 
discussed in Sec. Ill Al Hereafter, for reasons to be ex- 
plained shortly we neglect small correction it/(6ujht) <C 
1 in the screened electron-electron interaction 1121) . 



D. Hartree-Fock decoupling 

Effective action (|10|) involves electron states with spin- 
up and spin-down projections. Electron-electron inter- 
action can flip electron spin. Therefore, a charge den- 
sity wave state is characterized by an order parameter 
ApJp 2 2 (Q) that is matrix in the space of Landau level 
and spin indices. However, as it will be clear from dis- 
cussion below, if the Landau levels are spin-resolved, i.e. 
A cx ^ max{T, r^ 1 }, the charge density wave state cre- 
ates only on the Nth Landau level with certain spin pro- 
jection. Then Landau levels with different spin projec- 
tion become completely separated and can be ignored. 
Thus, we can consider the charge density wave order pa- 
rameter to be matrix only in the space of Landau level 
indices. It is related with distortion of electron density 
on the (N - l)th, Nth, and (N + l)th Landau levels as 

(p{q))=Sn L E A PiP2(q)^ P iP2(q)> (13) 

Pl ,P2=N—1 

where S stands for the area of two-dimensional electron 
gas and form-factor F P1P2 (q) is defined as 



^Pi P 2(q) = ™L 1 X^pifc(°) <7 Wq^)cx P ( ^qxq y i 2 H 



(14) 

After Hartree-Fock decoupling of interaction term in 
effective action (f 101) (see Refii&), we obtain 
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where 



dl "E E fe B ( r )[(««'n + A i 



N r n A 

T 

Pl ,P2 

-V dis (r))5 nm -H + \ PlP2 (r) C-J r )> ( 15 ) 



n A = 



n L S 2 



(16) 



Potential A PlP2 (r) in Eq. (|15|) appears as a consequence 
of distortion of uniform electron density by the charge 
density wave and is related with the order parameter as 

A PlP2 (q) ^E Tff vt")' 

pi P'? 1 qj 

P3P4 1 F 

where U PlP2P3Pi (q) denotes the generalized Hartree-Fock 
potential 

Up lP2 p 3Pi (q) = -n L (u scr (q)F plP2 (q)F P3Pi (-q) (18) 



dp 



(2i:) 2 n L 



- iqpl2H U SCT (p)F plP4 (p)F P3P2 (-p] 



E. Average over disorder 

After standard average over the random potential 
Vdis(r) (see Ref^), effective action (|TS)l becomes 

S = + f dr^(r) (ilu + n-n + X + iQ) ^(r) 



T 



~ J drtrQ 2 (r), 



(19) 



where we introduce new field Q(r), that is unitary matrix 
in Matsubara and replica spaces. For convenience we use 
the following notation 



= E E ^uw^w^w 



(20) 



Let us recall that action (|19|l at zero temperature, i.e. 
for ui n — > 0, and in the absence of the induced potential 
A(r) and the time-dependent vector-potential a has the 
following saddle-point solution 



(Ps P ) 



a/3 



P n S S' 



(21) 



where V is arbitrary global unitary rotation and P s ™ 
obeys the equation 



N+l 



P™ =igG^{v,v), G""(r,r') = E G p"( r < r ')- 

p=N-l 

(22) 

Here Green function Gp™(r,r') is as follows 



(23) 



Gp(u> n ) = [iuj n + MiV + £N ~ £p + iP, 



spj 



where chemical potential ^jv is measured from the A^th 
Landau level. The e p = uih(p +1/2) and p fc(r) are 
the eigenvalues and eigenfunctions of the hamiltonian 7i, 
and k denotes pseudomomcntum. In the case of weak 
disorder, ljjjt ^ L solution of Ea. (|22|l yields 



P 



sp 



signer, 
2r 



1 

2r 



'griL, 



(24) 
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with til = l/2irl 2 H . 

The fluctuations of the V field are responsible for the 
localization corrections to the conductivity (in the weak 
localization regime they correspond to the maximally 
crossed diagrams). However, in the considered case, these 
corrections are of the order of In N/N -C 1 and, therefore, 
can be neglected. For this reason we simply put V = 1. 

The presence of the induced potential A(r) and the 
time-dependent vector potential a results in a shift of the 
saddle-point value (|24|l due to the coupling to the fluctu- 
ations SP = P — P sp of the P field. The corresponding 
effective action for the SP field follows from Eq. H16|) 
after integrating out fermions: 



S = / rfrtrlnG" 1 - 



N r n A i 



T 2g 
drtrln 1 + (iSP + K + X)G 



drtr(P sp + 5P) 2 
(25) 



Finally, the thermodynamic potential can be written as 

fi = In J V[SP]I[6P] expS, (26) 

where following RefiSi the integration measure I[5P) is 
given by 



In I [SP] 



af3 



£[l-e(nmpP™«££,, (27) 



with p being the thermodynamic density of states and 
<d(x) the Heaviside step function. 

The quadratic in SP part of the action (|25|l together 
with the contribution 127|) from the integration measure 
determine the propagator of the SP fields (see Ref^ for 
details) 

_ gS mim J m2m3 S aS S^Q(mim 3 ) _ 2 [1 - 9(mim 3 )] 
1 + gn^i (w„ i3 - w mi , q) (irp) 2 
nS S a @ nS S 5 f 

X l + .g7r-™i(0,q) 1+ fl7 r^3 (0, q) ' 

where the bare polarization operator 7r" m (f„, q) involves 
Green functions for the (N - l)th, Nth, and (N + l)th 
Landau levels only 

7r Wm (v n , q) = Kr P 2 K> q) = - u l E G P2 ( w »0 

PlP2 P\P2 

x G Pl (aj m + Vn)Fp lP2 (q)Fp 2 p 1 (~q)- (29) 



F. Thermodynamic potential. Second order 
contribution 

In the absence of the time-dependent vector potential a 
effective action (|25H should contain only Ajvjv(q) = A(q) 
order parameter in the limit max{ T,t -1 } <C A cx <C u>h- 



To demonstrate it, we find the second order contribution 
to the thermodynamic potential for a = 0. 

Performing evaluation similar to one presented in 
Ref^, we obtain 



Q = si (0) + + . . . , 



where 



dr-trlnG -1 - 



1 



dr tr Pd 



(30) 



(31) 



2gJ sp 

is the thermodynamic potential of the liquid state and 



& 2 



n L S 2 
2T 



E 



dq 



T E E ^ 



P2P5P6 



i(q) 



(32) 



P3P4P7P8 



(-q) 



P5— P8 



F P5 p 6 (q) Fp7 P 8 ( q) 



x ^ P5P8 W- 1 + g ^„ (0;q) J^ 5P6 (0,q) 

x A PiP2(q) A p 3P 4(-q) 

is the contribution to the thermodynamic potential 
quadratic in the order parameter A plP2 (q). 

It is worthwhile to mention that polarization operators 
7r piP2 (^ n > q) °t> e y the following hierarchy with respect to 
small parameter, max{T, t^^/loh <C 1, 

'O(l), Pl=P2 = N, 

° ( ^f^ ) ■ Pi = ^orp 2 = iV, 

, ([ ' n "S T " 1} ] a ) 1 Pi^and^AT. 

(33) 

According to the hierarchy l|33|l we can write 

n PiP2 (°> q) ~ ^o" (°) <1)$pinS P2 n, (34) 

where we introduce 7To"(0, g) = 77^-^(0,(7). Thus from 
Eq.lO we obtain 



P1P2 



n L S 2 
2T 



E 



dq 

(2^ 



^ 1 + qtt; 



1---P4 

n L G 2 N (uj n ) 



+ .9<"(0,q) 
x A plP2 (q)A P3P4 (-q). 



^PlP2P3P4 

(q) (35) 

^pip 2 A r Af(q)^P3P4A r A r ( _ q) 



To find the possible non-zero order parameters 
A plP2 (q), we should diagonalize the 9x9 matrix in 
Ea. (|35|) . Fortunately, non-trivial part of Ea, 135|) can be 
written as 



a(q) a{q) 
a(q) a(q) + 2^(q) / 



A(-q) 



<^(-q) 



(36) 
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Here y(q) involves a linear combination of all order pa- 
rameters A pip2 (q) except A(q). Characteristic energies 
To(q) and Ti(q) is related with Hartrcc-Fock potential 
lfT5)l as follows 



T (q) 



UlVNNN(q) 



Ti(q)=e 



i<^t7jv,jv±i,jviv(q) 



(37) 

where denotes angle of vector q with respect to the x 
axis. We emphasize that quantity T\{q) depends only on 
the absolute value q of vector q. Matrix element a(q) is 
given by 



a(T, T -\q) = 1 + 4TV "^^M 

whereas function £(g) is defined as 

T Q (q) 1 



2T x (g) 



(38) 



(39) 



The eigenvalues of the 2x2 matrix in Eq. l|36|) can be 
easily found 



X± (q) = a (q) + £(g) ± VK^F+KWP 



(40) 



As one can check, the eigenvalue A+(g) has the same 
sign as £(g) for all values of a(q) whereas the eigenvalue 
A_ (q) changes its sign at point a(q) = 0. Therefore, the 
instability appears at the same condition as if we consider 
only one charge density wave order parameter A(q) as it 
has usually done According to the result derived 

in AppcndixlAl characteristic energy T\ (q) is of the order 
of T (q)/N < T (q). By using the condition £(g) > a(q), 
we find therefore 



TILS' 2 

2T 



dq 



Mq) 



a(q) 1 



Tx(q) 
Mq) 



a(q) 



A_(q)A_(-q) 



Tl(q) 
To(q) 



A + (q)A + (-q) 



(41) 



where 



A+(q) = ¥>(q) + o(?)A(q) 



a(gMq), (42) 



Minimum of the free energySi is reached at A+(g) = 0. 
Neglecting the difference of the order of 0(N~ 2 ) between 
A_(g) and A(g), we obtain finally 



ulS 2 



2T J (2tt) 2 
x A(q)A(-q). 



^ T (q)a(q)(l-^la(q) 



To(q) 



(43) 



Thus, the fact that the order parameters A pip2 (q) with 
Pi and P2 different from A can exist leads to correction 
of the order of 0(N^ 1 ). Later on we assume therefore 
that 



A PlP2 (q) = A(q)S PlN S P2N . 



(44) 



G. "Three-level" model 

The results of the previous section allows us to estab- 
lish finally an effective action for the "three-level" model. 

According to definition (|17|l . charge density wave on 
the Ath Landau level with the order parameter A(q) re- 
sults in the induced potential Ajv,jv±i(q), scattering elec- 
trons from the Ath Landau level to the (A± l)th Landau 
level. However, the induced potential is of the order of 
T\ (q) and, consequently leads to small corrections of the 
order of 0(N^ 1 ). For the reasons to be explained shortly, 
we write 

A PlP2 (q) = SU(q)F N 1 (q)A{q)S PlN S P2N . (45) 

Finally, the effective action for the "three-level" model 
becomes 



Stl[SP] = J drtrlnG 



T 



i J drtr(P sp +SP) 2 (46) 
drtrln[l + (iSP + K + P n XP n )g] , 



where 

Pjv(ri,r 2 ) = ^^ fe (r 2 )0ivfe(ri) 



= » £ e»p( >-^ 1+ ^ ) (47) 

- |ri :/ 2 n ^ nri ' r2 ' 2 



4/2 



II 2 



is the projection operator on the Ath Landau level 
(Ln(x) denotes the Laguerre polynomial) and 



n A = 



n L S 2 f dq 



J ^C7(g)A(q)A(-q), (48) 



Here, for a brevity we introduce U(q) = Unnnn{q)- 



III. CONDUCTIVITY OF THE UCDW STATE 

ATT c -r«T c 

A. Conductivity tensor a a b 

Effective action (|46|l allows us to evaluate conductivity 
of the system in the CDW state. As the most interesting 
case we consider the half-filled Ath Landau level where 
at T < T c the UCDW state exists. The order parameter 

:„18 



A(q) = 



S 



A 



*(q-Qo) + *(q-Qo) 



(49) 



where vector Q that determines period and direction 
of the UCDW state can be oriented along spontaneously 



6 




FIG. 1: Diagrams for corrections to the conductivity tensor 
dab- Solid line denote Green function, N/p/p' Landau level 
indices, dashed line the induced potential A(r), and shaded 
block the impurity ladder. 



chosen direction. Usually, its direction is fixed either by 
intrinsic anisotropy of the system or by small magnetic 
field applied parallel to 2DEG 22 i 23 i 24 i 25 i 26 . We assume 
that the vector Q is directed at an angle <j> with respect 
to the x axis. Let us recall that the absolute value of the 
vector Q equals Qo = ro/R c with tq « 2.4 being the 
first zero of the Bessel function Jo(x). 

The conductivity tensor <J a b(y n ,q) at q = can be 
found after integration over 8P{r) fields as the second 
derivative of logarithm of the effective action with re- 
spect to spatially constant time-dependent vector poten- 
tial a(V„), 



ttT 



SN r v n 8a a {v n )5ab{—v n ) 
In J V5PI[5P] expS TL [8P} 



(50) 



a=0 



It is worthwhile to mention that Eq. H50|) corresponds to 
the term j(z/ n )a(— v n ) with }(f n ) being current density in 
the effective action for the vector potential a(i/„). As one 
can check by inspection, contribution to the conductiv- 
ity tensor of the first order in the order parameter A(q) 
vanishes. It occurs because the UCDW state appears at 
non-zero vector Q . Thus, the first non- vanishing con- 
tribution to the conductivity tensor of the UCDW state 
is of the second order in the order parameter A. 

In order to find it, we expand the effective action 
Stl[5P] to the second order in both the induced po- 
tential A(r) and the K. Then, we integrate over <5P(r) 
fields. We do not present the explicit calculations here 
since they are similar to ones presented in the Ref>i£. 
We mention that there are three contributions of differ- 
ent structure to the conductivity tensor of the UCDW 
state. Diagrams for them are shown in Fig. ^ 

The first and the second diagrams (Fig. ^a)) corre- 



+ 



P, 



FIG. 2: Impurity ladder. Frequency u n + v n runs from right 
to left, whereas ui n runs from left to right. 



spond to the following contribution 

a « b[Vn) ~ ~^T J ° A 7 1 + ^(0,0) 

D a Np D b pN G p (uj n + v n ) 



(1+.9<"(0,Q )) 2 



(51) 



Here D^q p denotes matrix element of the covariant deriva- 
tive 



D a Np = I dr4,* Nk {r)D a <j) pk {v) = y/nZ^N-iFVN 

(52) 



where 



+ 5 p , N+1 -y a VN+ 1 



j x = i, 7* = 1, P x = -i, (3 y = 1. (53) 



We note that Eq. (|51J) contains only isotropic contribution 

to the conductivity tensor, and 

—Vyx for some direction of the UCDW. 
The third diagram (Fig. ^b)) is given by 

crffK) = ^ 2 A 2 r££- 



v n m ~ ' *—f (1 + g7To"(0, Qo)) 2 

x D pN D h Np ,Gp{u) n )G p i{u) n ) In pp 'n{Qo) 



x In p nn{Qo)Innp'n(Qq] 



(54) 



Symbol -f Pl p 2 p3P4 (Q) denotes the impurity ladder in the 
Landau level index representation (see Fig. |2J) 

/ PiP2P3 P 4(Q) = 9 1 F pip2 {q)F p3Pi (-Ci)cyLV (-iqQl 2 H ) . 

Jq 

, (55) 

Evaluation of Ip lP:2 p 3 p 4 (Qo) is given in Appendix iBl 
For convenience we present the results for quantities 
Jp iP2 p 3P4 (Q ) in the Table [fl As one can see, the 
Innn p n(Q ), In p nnn(Qo) and ^JVpp'w(Qo) arG propor- 
tional to J7o(t"o) = 0. Thus, the contribution (|54|l van- 
ishes 



=0. 



(56) 
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7jv,jv-i,jv-i,jv 


— 


gn L j£(ro) 


In,n+i,n+i,n 


= 


gn L Jo(ro) 


In,n-i,n,n+i 


= 


gn L Ji(r ) 


In,n+i,n,n-i 


= 


gn L Ji(r ) 


In,n,n-i,n 




giiLe"^ J (r )Ji{r ) 


In,n+i,n,n 




gn L e~ l,p Jo(ro) Ji(ro) 


In,n-i,n,n 




-gn^e* J (ro)Ji(ro) 


In,n,n+i,n 




-gnLe- i4, J (r )Ji(r ) 


In,n-i,n,n-i 




gn L e 2i *J?(r ) 


In,n+i,n,n+i 




gn L e- 2i *J?(r ) 


In,n-i,n+i,n 




gn L e 2i, l>Ji(ro)Jo(ro) 


In,n+i,n-i,n 




gnLe- 2i4, Ji(r )J (ro) 



TABLE I: Expressions for quantities I P1P2 P3P i (Qo) involved in Egs. l51H . l54H and I57H . 



The last diagram (Fig. ^c)) can be written as 



and l|59|) is seemed to be proportional to (T c — T)/T c 
since in the Landau theory the order parameter A oc 



n + Vn) (T c — T) JT~ c (see Ref^— for explicit expression) 

a ab \ Vn ) ~ 17~^7 2— t 2-^1 (1 +5 ^(0,Q )) - — - - 



G p (uj n )G p '(uj n + v n ) 
(l + 37r-»+^(0,Q )) 

D pN D p'N I NpNp' (Qq) 
1 + gKo"(v„, Qo) 



(57) 



We note that terms in the sum over Landau level in- 
dices with p = p = N ± 1 lead to the anisotropic contri- 
bution. Terms with p = N ± 1 and p — N =p 1 result in 
the isotropic contribution. 



B. Anisotropic contribution o^ ms ' 

We start our analysis of general expressions obtained 
in the previous section from Ea. l|57(l that contains the 
anisotropic contribution to the conductivity tensor. Tak- 
ing into account only terms with p = p = N ± 1 in the 
sum over Landau level indices involved in Eq. I|57|) we ob- 
tain 



(anis) 
<7 xx 
(anis) 



and 



(anis) 
&xy 



4nNJ?(r )h 



1 



4TTT r T 



Function h(z) is given as 



h(z) = 5z E 



C(M + z ) 
[C(2,i + ^)] : 




A 2 cos[2</>], (58) 



A 2 sin[2$. (59) 



z < 1, 



z > 1, 



(60) 

where £(fc, z) denotes the generalized Riemann zeta- 
function^l. Function h(z) increase monotonically from 
to 1, as it is shown in Fig. [3] 

Eas. (|58|) and l|59|l constitute one of the main results 
of the present paper. Anisotropic contributions l|58|) 



The (JxT as the function of angle <f> has the minimum 
at <fi — that corresponds to the vector Q directed along 
the x axis. We note that the modulation of electron 
density along the y axis is absent in this case. At <p = 
conductivity Oy%™^ is positive whereas a x a x nm ^ is negative, 
moreover, they have the same absolute values. It leads 
to the statement that the conductivity a xx (along the 
electron density modulation) is less than conductivity a yy 
(across the modulation). Wc emphasize that at <fi = 
conductivity a x a y ms ^ vanishes. 

If the vector Q is oriented at angle (j> = n /A with 



respect to the x axis conductivities a x a x ls ^ and a 



(anis) 



vanish due to the symmetry between the x and y axes. 
Vice versa, conductivity criy ms ^ reaches the maximum. It 
is worthwhile to mention that anisotropic contributions 
to the conductivity of UCDW state are proportional to 
N as the conductivity of the liquid statai. 



Isotropic contribution cr^ ct ' 



Eas. (|51|) and (|57|> allows us to find also the isotropic 
contribution to the conductivity tensor at T c — T <C T c . 
Taking into account Eq. (|51|l and terms with p = N ± 1 
and p = N =p 1 in the sum over Landau level indices in 
Eq. (I57|) , we obtain the following isotropic contribution to 
conductivity a xx 



1 



4irT r T 



(61) 



where function h xx (z) is given by 
hxx(z) = J 1 2 (r )h(z) + 



2zC(2, \ + z 
_ Rey/ (\ + (1 -i)z) 
C(2,i + z) 



(62) 
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FIG. 3: Functions h(z), h xx (z) and h xy (z) 



Here ifi(z) stands for the Euler di-gamma function. The 
h xx {z) has asymptotic expressions as 



3*' 



Ji(ro) + t - 



1 3J?(ro)-l 



z < 1, 



z > 1. 



(63) 

As one can see from Ea. (|61|l . isotropic contribution 

^cri'l *' to conductivity a xx is of the same order as the 
anisotropic contribution. 

The isotropic contribution to conductivity <j xy is as 
follows 



y uj h \AitT c t 



(64) 



where 



5 ' \ + z ) + Z ' C i^\ + z 



1 



2h !Eai (z) + zlmV (- + (l-«')2 



Im^> ( — + (1 — i)z 



[C (2, *+*)]' 



(65) 



The function h xy (z) has the following asymptotic expres- 
sions in the limits of small and large z 



h X y (z) 




Z > 1. 



(66) 

We mention that the isotropic contribution criy *' con- 
tains additional small factor max{T c , t -1 }/^^ compared 
to the others. Results <|61|) and (|64fl are one of the main 
results of the paper. 



IV. EFFECT OF THE ORDER PARAMETER 
FLUCTUATIONS ON THE CONDUCTIVITY 
TENSOR 

A. Order parameter fluctuations 

The order parameter A(r) has meaning of the saddle- 
point solution for a plasmon field that appears in the 
Hubbard-Stratonovich transformation^ of the screened 
electron-electron interaction. The expansion of such 
physical quantities as free energy and linear response 
in the order parameter series can be justified if fluc- 
tuations of the order parameter can be neglected. As 
it was showni^, fluctuations of the order parameter re- 
sults in the first order transition from the liquid state 
to the UCDW state at temperature T c — ST C where 
8T C /T C oc N~ 2 / z <C 1. In the present section we inves- 
tigate the effect of the fluctuations on the conductivity 
tensor above and below the mean-field transition. 

In the previous section we assumed that the direction 
of the CDW vector Q is fixed by intrinsic anisotropy of 
crystal or by applied parallel to 2DEG small magnetic 
field. However, the functional dependence of anisotropy 
term in the hamiltonian was insignificant for mean-field 
results obtained above. Now it should be concretized. 
Experimental research of the anisotropy that determines 
the direction along which the UCDW creates has been 
performed in a number of papers^^* 2 ^***^. The re- 
sults obtained can be explained if one suggests that the 
Hartree-Fock potential U (Q) involves terms proportional 
to cos 2(f) and cos A(f>. We mention that the term cos 2(f) 
can be derived when small magnetic field parallel to 
2DEG applied^. However, without parallel magnetic 
field the term cos2</> is restricted by the symmetry of 
bulk GaAs crystal. To date its physical origin is un- 
known^. As experimentally proven^, coefficient of the 
cos 2<f) term depends on the density n of electrons. More- 
over, at some certain value n* of the electron density it 
vanishes and next term proportional to cos 4</> becomes 
important. Below we restrict our discussion to the gen- 
eral case n ^ n*. We note that the typical value of the 
anisotropy energy Ea is of the order of 1 mK per electron 
as it is obtained from experiment^. In order to take into 
account the anisotropy quantitatively we perform the fol- 
lowing substitution (see Eq.ljSJJl) 



T (Q)^T (Q) + E A 



1 — cos 2<f> 



(67) 



near Q = Qo- We note that the expression above has 
minimum at = 0. 

At T > T c the UCDW order parameter is zero in av- 
erage (A) = but the average of its square is non-zero 
(A 2 ) 7^ 0. It results in additional contribution to the 
conductivity tensor of the liquid state. It is worthwhile 
to mention that the contribution discussed above is anal- 
ogous to one for normal metal due to superconducting 
paring above critical temperature^. 
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The additional contribution to the conductivity ten- 
sor due to the order parameter fluctuations can be found 
with a help of the substitution (A(Q)A(— Q)) for A 2 in 
Eas. H58ft . 15911 . and l|64|) and averaging over all possi- 
ble vectors Q. The Green function of the order parameter 
is as follows (see Ref^) 



(A(Q)A(-Q)) = 



4:T (Q )n L 



T-Tq 



1 



4TTT r T 



x (Q -QafRl+Vsin 2 



(68) 



where dimensionless parameter rj = Ea/Tq(Qo) and we 
introduce 




FIG. 4: Functions H(z), H xx (z) and H xy (z) 



(69) 



Here constant (3i « 2.58. After integration over absolute 
value of vector Q we find that in Eas. (|58|) . (|59|l . I|61[) and 
(|64|) the following substitution should be used 



We note that the function h(z) has monotonic growth 
whereas function £ (2, i + z) / \/j(z) monotonically de- 
creases. As a result the function H(z) has maximum at 
z ps 0.97 (see Fig. 0J. Function Fa(x) involves complete 
elliptic functions of the first and second kind 



2tt 



/ 1 



V 4ttT c t 



T-T r . 



f(4>) 



2tt 



VT C . 2 
■ sin - 



T-Tr. 



-1/2 



(70) 



F A {x) 



1 + - ) K(iVx) - ~E(iy/x) 

X J X 
X < 1, 

— = lnl6e _4 a;, x 3> 1. 

x 



(73) 



It is worthwhile to mention that in order to obtain the 
result for T < T c from the known result for T > T c we 
should substitute 2(T C - T c ) /T c for (T - T c ) /T c as usual. 



B. Fluctuational correction to the anisotropic 
conductivity cr^ m8 ' 

Integrating over angle <j) in Ea. (|70|) for /(</>) = cos 20, 
we obtain the following fluctuational corrections to the 
anisotropic part of the conductivity tensor above and be- 
low T r 



Integration over angle (f> in Eq.(7U|) with f(<fi) = 
sin 20 vanishes. Thus, fluctuational corrections to the 
anisotropic part of a xy and a yx conductivities are absent 
if the UCDW is oriented at the angle = with respect 
to the x axis. We mention that in this case the mean- 
field contribution to a xy and a yx vanishes as well. Conse- 
quently, the off-diagonal components of the conductivity 
tensor are isotropic for = 0. 

It is worthwhile to emphasize that Eq. (|71() constitutes 
one the main results of the present paper. 




1 



iirT c T 
r[T c 



T r . 



V2 *\2(T C -T) 
F A (^-' 

1 T-Tr 



T-T c 
, T < T c , 

T > Tr. 



(71) 



Here function H(z) is determined by the function h(z) 
as 



2tt 4 



H(z) 



((2,± + z)h(z) J 3v W 



V3 



z^J?{r )+Wi 



z < 1, 
z > 1. 



(72) 



C. Fluctuational correction to the anisotropic 
conductivity <J^ oti 



Integrating over angle in Ea. (|70fl with /(0) = 1, 
we obtain the following fluctuational corrections to the 
isotropic part of the conductivity tensor above and below 

T r . 



1 



4ttT c t 
VT< 



T r . 



T-T r 



(74) 



^2 1 \2(T C -T) 
F ( vTc ' 

1 \ T - T r . 



, T<T C1 



T>T r . 
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and 



g (isot-f) 



2nr T c ( JL 



f 1 



x < 



F, 



\4ttT c t 
rjT c 



T-T c 
, T<T C , 



(75) 



V2 \2(T C -T) 



Here functions H xx (z) and H xy {z) can be expressed via 
functions h xx (z) and h xy {z) as follows 



H xx (z) 



((2,± + z)h xx (z) 



----- 



H XV (Z) 



V3(l + 4,7?(ro)) l 



C(2,i + z)/i xjy (z) 

1 ^ 7^-2^(1) ^ 



2 < 1, 
2 > 1, 



(77) 



37T-8J-2M 
,4>/3^(ro) + 3/3i' 



2 < 1, 
z > 1. 



We mention that functions H xx {z) and H xy {z) have max- 
imum at z equal to 0.34 and 0.16, respectively, as it is 
shown in Fig.0J Function Fi(x) involves complete ellip- 
tic integral of the first kind 




0.20 -0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 0.20 

Temperature, (T-T ) / T 



(76) 

FIG. 5: The dependence of anisotropic part of conductivity 
{(jyy — a xx )/2a x °J on temperature for 1/4tcT c t = 0.24, r\ = 
0.01 and N — 2. Dashed line corresponds to mean-field result 



the facts that maxF/(x) = 1 and max. H xx (z) 
Fig. and Eq.JTSJ), wc obtain 



1 > 



\T C - T\ 



0.3 (see 



(80) 



As one can see, Ea. (l80ll is fulfilled in the limit N ^> 1. 
Analysis of Ea. (|71|l results in similar non-equality with 
0.1 instead of 1 in the right hand side. Ea. (|76|l contains 
additional small factor T c /loh- Therefore, the condition 
of applicability for results (|71f) . I|75|) and (|76|l is given by 
Eq.Q. 



V. DISCUSSION 



F J( ,) = ^=M hlfa] HI (78) 



We emphasize that Eas. (|76|) - (|77ll are one of the main 
results of the present paper. 



D. Limit of applicability of Eas. ffl} . lT75l and f77| 

Eas.l|7T]). lf73|) and have singularity at T — > T c . 
It indicates that the results are not applicable near T c . 
The limit of their applicability is determined by the con- 
dition that fluctuational corrections (JJTJ, (|75|) and ij7rj|) 
are still small as compared to conductivity of the liquid 
state equal toi 



r(°) 



2N 

7 

7T 



4°J = n. 



(79) 



Below we prove that the condition of smallncss of fluc- 
tuational corrections and the condition \T C — T\ <C T c 
are compatible. Let us first consider Ea. (|75|l . By using 



In the previous sections we derived a number of results 
for the conductivity tensor of two-dimensional interacting 
electrons on a half-filled high Landau level. Wc demon- 
strated that below temperature T c of transition from the 
liquid to the UCDW state the anisotropic part of the 
conductivity tensor emerges. At (T c — T)/T c <C 1 the 
anisotropic part is proportional to deviation of temper- 
ature from T c . As it is shown in Fig. [31 it results in a 
cusp of temperature dependence of the conductivity at 
T = T c . Fluctuations of the order parameter above and 
below transition temperature T c smooth out the cusp (as 
it is shown in Fig.[SJ). 

Results discussed above have been derived in the case 
of the white-noise random potential. In the case of a ran- 
dom potential with correlation length d arbitrary related 
with the magnetic length In, we can state that temper- 
ature dependence of the results above remain the same 
whereas the functions h a b and H a b become to depend not 
only on l/47rT c r but on ratio d/ln also. 

In experiments^ magnetoresistance R xx and R yy as 
functions of temperature at filling factor v = 9/2 have 
been investigated. Unfortunately, the detailed tempera- 
ture dependence near the point at which R xx and R yy 
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become to deviate from each other did not investigated. 
Nevertheless, results reported in Refsi^ confirm the lin- 
ear dependence of magnctorcsistance R xx and R yy on 
temperature in the certain range of temperatures not 
too close to the T c . Therefore, the detailed investiga- 
tion of temperature dependence of magnetoresistance are 
needed in the future. 

The results <|58|) and l|59|l for the angle dependence of 
the anisotropic part of the conductivity tensor qualita- 
tively describe the results of experiments&iiii. If a cur- 
rent runs in the direction of charge density modulation, 
the conductivity (resistance) is less (more) than the con- 
ductivity (resistance) for the case of current running in 
the direction perpendicular to charge density modula- 
tion. The Hall conductivity (resistance) remains roughly 
the same. 

The UCDW state exists also near half-filled high Lan- 
dau level. In the case of /ijv 7^ 0, the functions h a b and 
H ao & re dependent not only on parameter 1 /AttT c t but on 
ratio /hjv /T c also. Increasing the chemical potential [In 
we transform the UCDW state to the triangular CDW 
state characterized by vectors Q 1; Q 2 and Q 3 directed 
at angles <j>, <j> + 2ir/3 and <f> + 4ir/3 with respect to the 
x axis, respectively. However, due to the identity 



1 + e 2m/3 + e - 27rl/3 = 



(81) 



the contribution from three vectors to the anisotropic 
part of conductivity tensor vanishes. Thus, the con- 
ductivity tensor of the triangular CDW state remains 
isotropic in the approximation which is of the second or- 
der in the CDW order parameter. From general physical 
arguments it is clear that the conductivity tensor of the 
triangular CDW state should be isotropic. 

It is worthwhile to mention that from physical point 
of view the anisotropic part c^, nis ' ) of conductivity ap- 
pears due to the existence of the induced anisotropic 
potential A(r) in action I0SJ. Anisotropic resistance of 
two-dimensional electrons in a weak magnetic field in 
the presence of unidirectional periodic potential has been 
measured at several Kelvins in heterostructures with mo- 
bility fiQ ~ 10 6 cm 2 /Vs by Weiss, von Klitzing, Ploog, 
and Weimann fifteen years ago^I. Theoretically, the ef- 
fect of unidirectional periodic potential on the conductiv- 
ity tensor of two-dimensional electrons in a weak mag- 
netic field has been investigated with the help of both 
the kinetic equation 3 ^ 3 , and the diagrammatic tech- 
m q UQ 34 i 35 i 36 However, the case considered in the present 
paper did not analyze in the Refsi 32 i 33 i 34 i 35 i 36 . The im- 
portant difference of the induced potential A(r) from ex- 
ternal periodic potential is that the induced potential 
A(r) appears on the iVth Landau level only whereas ex- 
ternal potential scatters electrons on all Landau levels. 



VI. CONCLUSION 

We obtained the conductivity tensor of two- 
dimensional electrons in the presence of weak disorder 



and weak magnetic field at half-filled high Landau level 
where the UCDW state exists. In the framework of the 
order parameter expansion we derived that at T c — T <C 
T c anisotropic part of the conductivity tensor propor- 
tional to (T c — T)/T c emerges. Also we demonstrated 
that the order parameter fluctuations result in additional 
anisotropic contributions near T c to the conductivity ten- 
sor that wash out the mean- field cusp at T = T c . The 
results obtained are in agreement with the experimental 
findings. 
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APPENDIX A: CALCULATION OF 
CHARACTERISTIC ENERGY Ti 

According to Eas. (|18|) and (|37|l the T\(q) is given by 
Ti(-z) = -n L e l<t, [U scl .(q)MqRc)Ji(qRc)e^ (Al) 



dp 



>Px - Wy 



i^fn L e Usc ^ P) p 2 l H 

x L NN I — — I Ljv,jv— 1 I — y~ 

where we use the following result- 



F, 



N,N 



_x(q) = V2N 



1x ~ Hy -q 2 l 2 „/4 



L 



N,N-1 



n 2 ! 2 



q 2 l H ' "'" 'I 2 

« e'^MqRc), qR c ^2N. (A2) 

The characteristic energy T\ = Ti(Qq) is given by 



4iV 



TsU H f dx MANx)MANx^ 
T i = "T^T I n Ji(2r x), (A3) 



4v^2 J e{x) 




where 



xV2 y 



e(x) = 1 + -^={1 - J 2 (ANx)). 



(A4) 



Performing calculation of the integral, we find 



Ti 



1 



2 V s/2r r s J l + V2r Q r s 



ci 



(A5) 
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where constant c\ equals 
1 



dx 



7T J xJx{\ - X 2 ) 
l/2r 



3in \ 2tqx 



097. 



(A6) 



As we can see from Ea. l|A5(l the characteristic energy 
T\ ~ Tq/N as we mentioned above. 



APPENDIX B: CALCULATION OF THE 



l(Qo) 



we find 



In,n-i,N-i,n = In.n+i,n+i.n 

2 f, Jo(2r x) 
= gn L - dx 

7T J VI -X 2 



l N,N-l,N,N+l — J iV,jV+l,/V,/V-l 

1 



(B3) 



2 f 

griL— / tix(l — 2x ) 



Jb(2r a:) 



Using definition (|55|1 and Ea. i|A2(l . we obtain in the 
limit N > 1 



IN,N,N,N-1 



= -I 



N,N,N,N+1 



IN,N-1,N-1,N 



dx 

iN,N+i,N+i,N = gn L N I — e~ 

x 



7T J \Jl- X 2 



(Bl) 



In,n-i,n,n+i — In,n+i,n,n-i — 9 n L^ / — e 

X 



x Ljy (a;)L Ar+1 (a;)J (r v / 2x), 



In,n,n,n-i — —In,n,n,n+i = griL^yN 



dx 



I 



N,N+1,N,N 



— 1 N,N+1,N,N 

gn^e -1 ^ — [ dxx^ 1 

7T J 



VT^x 2 



x L- N \x)L N {x)J 1 {r V2i), 



LN,N+1,N,N 



= -I 



n,n+i,n,n — gn-Le 



dx 



x e x L N L +1 (x)L N (x)J 1 (r V2x) 



'N,N-1,N,N-1 



gnrJ** I ' dx J ^ X 



7T j yi - X 2 





7jv,jv-i,jv,jv-i — ^Iv.jv+i.jv.jv+i = gri^e^N / 



(7.r 
.r 



In,n-i,n+x,n = In,n+i,n-i,n ~ gnLe 2l ^N / — e 



With a help of asymptotic expression for Laguerre poly- 
nomial^ 



L%(x) 



1 



irx \x 



nN (2a-l)/4 



cos 2Jnx - 

2 4 



AT>1, 



(B2) 



l N,N-l,N+l,N 



T* 

1 N,N+1,N-1,N 

gn L e^ 2 -[dx(l-2x*)?p4 

IT J X 2 



The integrals can be evaluated by using the following 
equalitj^l 



tt/2 

^ /" d4>cos(2fi(f))J2u(2rocos(t>) = J v +^(ro)J v -^(ro). 
o 

(B4) 

Finally, it yields the results presented in Tabic [I] 
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